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ABSTRACT 

The nature of the X-ray source RX J1914+24 has been the subject of much debate. 
It shows a prominent period of 569 sec in X-rays and the optical/infra-red: in most 
models this has been interpreted as the binary orbital period. We present our analysis 
of new XMM-Newton and Chandra data. We find a longer term trend in the XMM- 
Newton data and power at 556 and 585 sec in 5 sets of data. It is not clear if they 
are produced as a result of a beat between a longer intrinsic period and the 569 
sec modulation or if they are due to secular variations. We obtain a good fit to the 
XMM-Newton spectrum with a low temperature thermal plasma model with an edge 
at 0.83keV. This model implies an unabsorbed bolomctric X-ray luminosity of 1 x 
10 33 ergs s _1 (for a distance of lkpc) - this is 2 orders of magnitude lower than our 
previous estimate (derived using a different model). If the distance is much less, as 
the absorption derived from the X-ray fits suggest, then it is even lower at ~ 3 x 10 31 
ergs s -1 . 
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1 INTRODUCTION 

The nature of the X-ray source RX J1914+24 has been the 
subject of much debate in recent years. Only one modula- 
tion period (569 sec) has been seen in both its X-ray and 
optical light curves (Ramsay et al 2000). For this, and other 
reasons, the period of 569 sec has been taken to be the bi- 
nary orbital period, which would make it the second most 
compact binary system known. 

The models which have been put forward to explain the 
observational characteristics of RX J1914+24 fall into two 
broad categories: accretion and non-accretion driven mod- 
els. The fact that the optical spectrum of RX J1914+24 
(also known as V407 Vul) shows no emission lines, puts a 
question mark over the accretion driven models. The non- 
accreting model is the unipolar-inductor (UI) model of Wu 
et al (2002). While this model has been able to explain many 
of the observed characteristics of this source, it appears that 
its X-ray luminosity implies an induction torque which is 
much greater than the amount generated by gravitational 
radiation. There is some evidence that the observed spin-up 
rate of RX J1914+24 (Strohmayer 2002, 2004b, Ramsay et 
al 2005a) is therefore much smaller than that predicted by 



the UI model (Marsh & Nelemans 2005). This conclusion 
is based on the X-ray luminosity determined using ROSAT 
(Cropper et al 1998) and XMM-Newton observations (Ram- 
say et al 2005a) . The ROSAT data were well fitted using an 
absorbed blackbody model while the higher signal to noise 
data of XMM-Newton found a poor fit using this model. 
Ramsay et al (2005a) found the best fit to the data could be 
achieved using an absorbed blackbody plus Gaussian model, 
although they did not expect that this model was physically 
realistic. 

Since those first XMM-Newton observations were made 
in Nov 2003, a further two sets of observations made in 
Sept/Oct 2004 have become publically available. These ob- 
servations were longer than the initial XMM-Newton obser- 
vations and the particle background was lower. In addition, 
a further set of X-ray observations of RX J1914+24 made us- 
ing Chandra has become publically available. We have there- 
fore extracted these data from the archive and use them to 
make a more detailed investigation of its X-ray spectrum, 
to determine if the 569 sec is continuing to spin-up at the 
same rate and to search for evidence for periods other than 
the 569 sec period. 
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XMM 


Start Date 


Duration 


Mean EPIC 


Orbit 




(ksec) 


pn (Ct/s) 


0880 


2004-09-28 


17.0 


1.08 


0882 


2004-10-02 


19.0 


1.16 


Chandra 


Start Date 


Duration 


Mean Chandra 


Obs Id 




(ksec) 


Ct/s 


4519 


2004-04-05 


15.4 


0.44 


4520 


2004-06-13 


14.7 


0.37 


4521 


2004-09-05 


15.2 


0.34 



Table 1. The log of the observations of RX J1914+24 reported 
here. Each observation was continuous with no data gaps. Top 
panel: The XMM-Newton observations of RX J1914+24 made 
in 2004. The earlier observations are detailed in Ramsay et al 
(2005a). Bottom panel: The Chandra made in 2004. The earlier 
observations are detailed in Strohmayer (2004b). 

2 OBSERVATIONS AND DATA REDUCTION 

2.1 XMM-Newton observations 

RX J1914+24 was observed twice using XMM-Newton in 
Nov 2003. The details of these observations has been re- 
ported by Ramsay et al (2005a). It was observed again late 
in 2004 and we show in Table 1 the start time, duration and 
mean X-ray count rate of the second epoch XMM-Newton 
observations. In this paper we discuss the data collected us- 
ing the EPIC X-ray detectors (pn and MOS) which are sen- 
sitive over the energy band 0.15-10keV. We extracted the 
data from the public archive and processed the data using 
the XMM-Newton SAS v6.5 in a similar manner to that de- 
scribed in Ramsay et al (2005a). Data were barycentrically 
corrected and in units of Terrestrial Time (TT). 

2.2 Chandra observations 

Observations of RX J1914+24 made using Chandra on 5 
occasions were reported by Strohmayer (2004b). Since then 
another 3 observations have entered the archive. All of these 
observations were made in continuous clocking mode: in this 
mode a one-dimensional image of the sky is produced. We 
have extracted all of the available data from the archive and 
corrected and analysed them in a similar manner to that 
carried out by Strohmayer (2004b). 



3 DETERMINING THE RATE OF SPIN-UP 

Ramsay et al (2005a) folded all available X-ray data 
(ROSAT, ASCA, Chandra and XMM-Newton) on the best 
fit initial period of Strohmayer (2004b) (ie without the pe- 
riod derivative) and found that the sharp rise of the X-ray 
flux arrives earlier in phase over time. A best fit to the 
data showed that RX J1914+24 is spinning up at a rate 
of 3.2xl0~ 12 s/s, which was similar to that reported by 
Strohmayer (2002, 2004b). We did the same analysis as done 
in Ramsay et al (2005a) and show the phase of the rise in 
X-ray flux for the Chandra data taken in April 2004 and the 
new XMM-Newton data in Figure 1. We find a best fit of 
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Figure 1. The phase of the sharp rise to maximum X-ray flux 
determined from ROSAT, Chandra and XMM-Newton data (see 
Ramsay et al 2005a for details of these observations). The more 
recent point refers to the new XMM-Newton data presented here, 
with the preceding point refers to Chandra data taken in April 
2004. 

P = -3.31±0.09 x 10~ 12 s/s (the error being determined us- 
ing a bootstrap method), which is consistent with the value 
we determined previously. 



4 SEARCHING FOR PERIODS OTHER THAN 
THE 569 SEC PERIOD 

Strohmayer (2004b) detected sideband structure in the 
power spectrum of RX J1914+24 in observations made us- 
ing Chandra but concluded that this was likely to have been 
caused by the dithering pattern of those observations. On 
the other hand, longer term variations have been seen in 
X-ray data (Ramsay et al 2000, Strohmayer 2004b). 

We have searched for evidence of periods other than 
the 569 sec period in these new X-ray data. We performed 
a Discrete Fourier Transform (DFT) on the combined data 
taken in XMM-Newton orbits 0718 and 0721 and also orbits 
0880 and 0882 as well as each of the Chandra observations. 
As expected, the 569 sec period and its harmonics were seen 
strongly in all the amplitude spectra. To pre-whiten the data 
we used two methods: we pre-whitened in the usual man- 
ner using a sinusoid at predicted period at the observation 
epoch (based on the initial period of Strohmayer (2004b) 
and the spin-up rate reported here). We also pre-whitened 
the data by subtracting the mean X-ray pulse profile ap- 
propriate to each individual observation. The latter is more 
appropriate since the X-ray pulse profile is non-sinusoidal. 
As expected most of the peaks are removed after this pro- 
cedure. However, in the second epoch set of XMM-Newton 
data (0880 and 0882), side-band structure is clearly present 
at 584.1 sec and 552.8 sec in the pre-whitened spectrum 
(Figure 2). They are located in the wings of the main peak 
in the original amplitude spectrum and have an amplitude of 
~1/10 of the main peak. In the earlier XMM-Newton data 
(0718 and 0721) there is only weak evidence of these peri- 
ods. In the Chandra dataset of 2004 June 13 we also find 
strong evidence for periods very similar to those seen in the 
XMM-Newton data (Figure 2). In two others, (2004 Jan 31 
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Figure 2. The pre- whitened amplitude spectra of 5 sets of X- 
ray data. They have been pre-whitened by subtracting the mean 
X-ray profile appropriate to each epoch. 



and 2004 Sept 5) there are peaks in the amplitude spectra 
which correspond to these periods, although they are only 
marginally significant compared to other peaks in the spec- 
tra. There is no evidence for these periods in the longest 
Chandra dataset which lasted 35 ksec. We also determined 
how the X-ray pulse profile changed over the individual ob- 
servation duration. We find that the mean profile brightens 
over the duration of both the 0880 and the 0882 observa- 
tions. 

We believe that 553 and 584 sec periods could result 
when a modulation signal (ie the 569 sec period) is mixed 
with another longer term variation (either a periodic or non- 
periodic modulation) . If it was a periodic modulation, then 
its period would be the difference between the 569 sec period 
and the 553 and 584 sec periods: just over 6 hrs. Alterna- 
tively if a modulation signal is mixed with either a 553 sec 
or 584 sec period then this would give rise to the longer term 
trend in the data. The question remains if this longer term 
trend is periodic or due to the secular variability that have 
already been observed in this system: we discuss this further 
in 86.1. 




channel energy (keV) 

Figure 3. The XMM-Newton EPIC pn spectra from the data 
taken in 2004. The solid line is the best fit to the model consisting 
of an absorbed thermal plasma model with an edge. 



5 SPECTRA 

Ramsay et al (2005a) showed the X-ray spectrum of RX 
J1914+24 was poorly modelled using an absorbed black- 
body model. Their best fit was obtained using an absorbed 
blackbody plus a broad Gaussian line near 0.6keV. This fit 
was still poor (x 2 =1.63) and the model was unlikely to be 
physically realistic. 

With the better quality data afforded by the new XMM- 
Newton data we have re-visited the X-ray spectrum of RX 
J1914+24. We fitted the EPIC-pn data from XMM-Newton 
orbits 0880 and 0882 separately. Given the unphysical nature 
of the previous fits and working from the satisfactory fits to 
our X-ray observations of other compact binaries (the AM 
CVn stars, Ramsay et al 2005b) we tried fitting the spectra 
with an absorbed thermal plasma model with variable ele- 
mental abundance (vmekal in XSPEC). This model still gave 
a poor fit to the data (x 2 =4.6 for the 0880 data). However, 
when we added an edge to the model we obtained a good 
fit to the 0880 data (x 2 =0.97, 122 dof; Figure 3). The fit 
to the 0882 data using the same model was slightly poorer 
(x£=1.21, 141 dof). Elements which were required in the 
fit include N, O, Ne, Ca, Fe for the 0880 data, with some 
evidence for Na in the 0882 data. The edge is located at 
0.83keV which is consistent with it being due to OVIII. The 
absorption derived from the fits is ~ 1.2 x 10 21 cm -2 . We 
note we cannot obtain a good fit with this model when we 
fix the absorption at 8.3 xlO 21 cm~ 2 , the larger value deter- 
mined from optical interstellar absorption features (Steeghs 
et al 2006). 

The observed flux in the 0.1-10keV energy band 
(1.30xl0~ 12 



ergs s 
1.48 xl0~ 12 ergs s _1 cm 



cm 

2 



2 for the orbit 0880 data and 
for the orbit 0882 data) is very 



similar to that found using the blackbody plus Gaussian line 
of Ramsay et al (2005a). However, the unabsorbed bolomet- 
ric flux in the fit presented here (1.14x 10 -11 ergs s _1 cm~ 2 ; 
0880 and 1.24xlO _n ergs s" 1 cm" 2 ; 0882) is substantially 
reduced compared to the earlier results (3.2 xlO -10 ergs s _1 
cm -2 ) where a blackbody spectrum was used. 
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6 DISCUSSION 

6.1 The periods in RX J1914+24 

Until now there has been no convincing evidence for any 
other periods in RX J1914+24 apart from the 569 sec period. 
In the case of 5 X-ray observations we have found evidence 
for power at 552 and 584 sec. This changes our perception 
of RX J1914+24. The data are clear - the X-ray flux clearly 
increases in the XMM-Newton data taken in both orbit 0880 
and in orbit 0882. This could be considered to be the result 
of an underlying ^6 hour period, or it could be the result 
of secular variability known already in this system at longer 
timescales (Ramsay et al 2000, Strohmayer 2004b). Both 
could produce the observed amplitude spectrum. 

To explore this question, we have obtained the mean 
X-ray flux in each of the Chandra and XMM-Newton ob- 
servations. Since the Chandra ACIS-I detector has a low 
effective area at energies less than ~0.4keV, we determined 
the observed mean flux in the 0.4-1.0keV energy band. We 
show how the flux in this band varies over time in Figure 4. 
We also show those epochs at which we detected the 552 sec 
and 584 sec periods which could be signatures of a longer 
period trend in the observation. Although not conclusive, it 
appears that at the higher flux levels there is no (or weak) 
evidence of these other periods, so that shorter (<1 day) 
timescale variations may take place only when the system 
is emitting at intermediate or low levels. The detection of 
power at 552 and 584 sec does not appear to be related 
to the observation length since the longest observation (the 
Chandra observation of Feb 2003 which lasted 35 ksec) did 
not show obvious power at these periods. 

At this stage, we cannot distinguish between the X-ray 
long term behaviour being due to random changes in either 
the accretion rate (in the case of the accreting models) or 
the electrical current dissipation (in the case of the unipolar- 
inductor model) . We urge the search for radial velocity vari- 
ations on a timescale of less than 1 day, or dedicated X-ray 
and/or optical observations of this system to search for vari- 
ations on this period. 

If the longer term variation in the X-ray flux is periodic, 
then we showed earlier that its period would approximately 
be 6 hrs. Smith & Dhillon (1998) show that for a binary 
system with this orbital period then its spectral type would 
be expected to be a mid K-type, although the range is up to 5 
sub-classes. We note with interest the discovery of Steeghs et 
al (2006) of features in the optical spectrum of RX J1914+24 
which appear to resemble a G9/K0 spectral type. Therefore, 
an orbital period of 6 hrs would be marginally consistent 
with this. 
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Figure 4. The observed flux in the 0.4-lkeV energy band (units of 
10~ 12 ergs s~ 1 cm~ 2 ) determined in Chandra and XMM-Newton 
observations. 'XP' denotes the presence of periods at 552 and 
583 sec: in the case of the most recent Chandra observation the 
significance is marginal - 'WXP' weak X-ray periods. 



the new model is an optically thin plasma, a correction ef- 
fect for projection is now no longer applicable. Moreover, 
absorption derived from the X-ray fits is ~15% of that de- 
rived from the optical data of Steeghs et al (2006) based 
on a G9 spectral type. If the G star is not the secondary, 
then we can reduce the distance by a similar factor and the 
luminosity is therefore 3xl0 31 ergs s _1 . In this case, the X- 
ray luminosity of RX J1914+24 is more than ~3 orders of 
magnitude less than that found by Ramsay et al (2005a). 

In the unipolar-inductor model of Wu et al (2002) the 
amount of power which is emitted is determined by vari- 
ous factors, with the degree of asynchronisity being one of 
the key determinants (along with component masses and 
magnetic field strength). However, to produce the observed 
luminosity, the asynchronisity would have to be 0.001 or less 
(Wu et al 2002). This reduces the objections to the UI model 
on energetic grounds raised by Marsh & Nelemans (2005) to 
a significant degree. 

The abundances are clearly non-solar, which is consis- 
tent for an evolved secondary as is seen in the fits to spec- 
tra of AM CVn systems (eg Marsh, Home & Rosen 1991, 
Strohmayer 2004a, Ramsay et al 2005b). If this model re- 
flects the properties of the material leaving the secondary 
star then it rules out the G9/K0 star seen in the optical 
spectrum being the secondary star. 



6.2 The X-ray spectrum of RX J1914+24 

We have shown that the derived luminosity of RX J1914+24 
is highly model dependent. Further, the distance to RX 
J1914+24 is still rather uncertain, but Steeghs et al (2006) 
suggest a distance of lkpc based on its optical spectrum. 
For a distance of lkpc the mean bolometric luminosity from 
the XMM-Newton orbits 0880 and 0882 data is therefore 
1.4x 10 33 ergs s _1 - an order of magnitude lower than the pre- 
vious luminosity of 3.8 xlO 34 ergs s _1 . Ramsay et al (2005a) 
went further and applied a correction to account for projec- 
tion affects giving a luminosity of ~ 1 x 10 35 ergs s _1 . Since 
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